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We consider Yang-Mills theory with the U{1) gauge group on a non-commu- 
tative plane. Perturbatively it was observed that the invariance of this the- 
ory under area-preserving diffeomorphisms (APDs) breaks down to a rigid 
subgroup SL{2,R). Here we present explicit results for the APD symme- 
try breaking at finite gauge coupling and finite non-commutativity. They 
are based on lattice simulations and measurements of Wilson loops with the 
same area but with a variety of different shapes. Our results are consistent 
with the expected loss of invariance under APDs. Moreover, they strongly 
suggest that non-perturbatively the SL{2, R) symmetry does not persist ei- 
ther. 
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1 Introduction 



Invariance under area-preserving diffeomorphisms (APDs) [1] is a basic sym- 
metry of ordinary Yang-Mills theories in two dimensions. In particular it 
means that Wilson loop expectation values only depend on the oriented ar- 
eas singled out on the manifold. Thanks to this property the theory acquires 
an almost topological flavour [2] and, as a consequence, it can be solved an- 
alytically. Elegant group theoretic methods [3] lead to closed expressions for 
the partition function and a set of observables [4]. 

The invariance under APDs was initially believed to persist also in U (n) 
gauge theories defined on a non-commutative (NC) two-dimensional mani- 
fold. It was assumed to play a central role in the large gauge group — 
characteristic of gauge theories on NC spaces — which merges internal and 
space-time transformations. A detailed study of the non- commutative gauge- 
transformation algebra was performed in Ref. [5]. If APD symmetry holds, 
one might hope to be able to solve gauge theories also on a NC plane by 
generalising the powerful geometric procedures developed in the commutative 
space. 

This scenario was suggested by an intriguing observation for U (n) gauge 
theory on a NC torus. It can be related by Morita equivalence to its dual on 
a commutative torus [6], where the APD invariance is granted. The theory 
on the NC plane would then be reached by a suitable limit, and one could 
hope for the invariance to be preserved [7]. 

Wilson loop perturbative expansions in the coupling constant g and in 1/6* 
— 6 being the non-commutativity parameter — were performed on the NC 
plane. To the order first considered in Refs. [8,9], the results were consistent 
with APD invariance. 

Later on Ref. [10] extended those results to the next order, namely 0{6~'^) 
at 0(5'^). The outcome revealed different expectation values for a Wilson loop 
with the shape of a circle and a rectangle of the same area. This observation 
motivated the systematic investigation in Ref. [11], where Wilson loops in 
a wide class of contours were considered in the axial gauge. These results 
suggest that the APD symmetry breaks down to a residual subgroup of linear 
unimodular transformations, SL{2,R). 

Subsequently non-perturbative arguments for this APD symmetry break- 
ing were given based on the Morita duality on tori [12]. Recently, Ref. [13] 
reconsidered this issue by applying twist deformation techniques, which also 
confirm that the APD symmetry may break at the quantum level. A new 
study of the large 6 expansion [14] reports the break-down of the area law at 
0(1/^^), in agreement with simulations results at finite 6 [15]. Ref. [14] does 
not report the observation of any symmetry. However, the question about 
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the ultimate status of (partial) APD symmetry is still open. 

The present work presents explicit results for Wilson loops with polygonal 
contours, at finite 9 and g, under APDs. Our non-perturbative results are 
obtained on the lattice and extrapolated to the continuum. They agree with 
the breaking of this symmetry, both on the lattice and in the continuum limit. 
Moreover we provide evidence against the survival of a residual symmetry 
subgroup SL{2,R). 

In Section 2 we briefly review the U{1) gauge theory on a NC plane, its 
lattice discretisation and the mapping onto a twisted Eguchi-Kawai model, 
which can be simulated. Section 3 presents our simulation results for the 
planar limit, which is necessary to identify a physical scale. Then we address 
in Section 4 the Double Scaling Limit to a continuous NC plane of infinite ex- 
tent, which allows us to study explicitly the effect of APD transformations on 
Wilson loops. In Section 5 we confirm the APD symmetry breaking by con- 
sidering observables, which differ from those in Section 4. Their perturbative 
treatment is commented on in an Appendix. Section 6 focuses specifically 
on the SL{2, R) symmetry, and Section 7 is dedicated to our conclusions. 

2 ^(1) gauge theory on a non- commutative 
plane 

In this work, we consider the simplest version of a Euclidean NC plane by 
assuming a constant non-commutativity parameter 6, so that the coordinate 
operators fulfil 

[x^,x^] = iOef,^ (/i,z/ = 1,2). (1) 

Such coordinates describe a charged particle moving in a (commutative) 
plane, which is crossed by a strong, orthogonal magnetic field. The lat- 
ter can be formally interpreted as i? oc 1/6, see e.g. Ref. [16]. A similar 
concept is also used to map open strings in a magnetic background onto NC 
field theory [17]. 

We can return to the use of ordinary (commutative) coordinates if all the 
fields are multiplied by star products (or Moyal products), 

(l){x) -k ^l){x) := 0(x)exp Q'9^6'e^^^j ^Ij{x) . (2) 

Here we focus on pure U{1) gauge theory with the Euclidean action 

S[A] = d^xF^^^F^, , 

F^,u = d^A^ - d^A^ + ig[A^, A^]^ . (3) 
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The last term is a star-commutator, which shows that even the U{1) gauge 
field is self-interacting on NC spaces. This action is star-gauge invariant, i.e. 
invariant under transformations 

A^{x) ^ U{x) * A^{x) * U{xy - ^U{x) * d^U{xy , (4) 

if U{x) is star-unitary, U{xy -kU{x) = 1. 

Other U (n) gauge theories may be studied along the same lines, but the 
formulation of SU (n) gauge theories runs into trouble on NC spaces. There- 
fore it is motivated to concentrate on f/(l) as a physical gauge group, which 
can be accommodated on NC manifolds. 



Although the points in such spaces are somewhat fuzzy, it is possible to 
introduce a lattice structure^ This is a first step towards a formulation to 
be used in Monte Carlo simulations. In the operator formalism this step 
imposes the constraint 

exp (i-^^m) = ^ y (5) 

where a is the lattice spacing. If we require the momentum components to 
be commutative and periodic over the Brillouin zone, the above condition 
implies that only discrete momenta occur, which is characteristic for a finite 
volume. On a. N x N lattice with periodic boundary conditions, the allowed 
momenta are spaced by 27r/ (aN). As a consequence, the non-commutativity 
parameter can be identified as 

e = -Na^ . (6) 

We are most interested in a Double Scaling Limit (DSL) 

a — > and N oo at Na^ = const. , (7) 

which leads to a continuous NC plane of infinite extent. The requirement to 
take the UV and IR limits simultaneously in a balanced way is related to the 
generic UV/IR mixing of the divergences in NC field theory [20]. 

This is clearly distinct from the planar limit, N ^ oo at fixed gauge 
coupling, which means here a fixed lattice spacing. The non-commutativity 
parameter diverges in this limit. In higher dimensions, this implies that non- 
planar contributions are suppressed, and the planar limit restores commuta- 
tivity in perturbation theory|§ In two dimensions the situation is different, 

"'^Here we only sketch this regularisation very briefly, for details we refer for instance to 
the review [18], or the theses quoted in Refs. [15] and [19]. 

^We remark, however, that this restoration does not need to hold generally: it can fail 
non-perturbatively for instance in the case of spontaneous symmetry breaking [21]. 
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and non-planar diagrams provide "anomalous" perturbative contributions in 
the limit of infinite non-commutativity, which are of the same order of mag- 
nitude as the planar diagrams [8-10]. However, such terms are shown to 
disappear when applying the procedure recently introduced in Ref. [14]. 

Even on the lattice it is far from obvious how to simulate NC gauge 
theory; note that the compact formulation seems to require star-unitary link 
variables. In this respect, it is highly profitable to map the system onto a 
twisted Eguchi-Kawai model (TEK model). This model is defined on a single 
space-point and its action takes the form [22] 

Stek[U] = -N(3 Z>.uTt (u^UMlUl'^ ■ (8) 

Ui and U2 are unitary N x N matrices which encode the degrees of freedom 
of the f/(l) lattice gauge theory. For the twist factor we adopt the choice of 
Ref. [15], Z21 = = exp{m{N + 1)/N) , where has to be odd. There is 
an exact equivalence to the lattice NC U{1) gauge theory, i.e. the algebras 
are identical, as Ref. [23] showed in the large limit. A refined consideration 
found such a mapping even at finite A^ [24]. Hence the TEK model can be 
used for numerical simulations of NC gauge theories, and it is most suitable 
for this purpose. 

It is straightforward to formulate Wilson loops in this matrix model. For 
instance, a rectangular loop with side lengths al and aJ (and clockwise 
orientation) corresponds to the term 

W{I X J) = Tr (ulU^Uruy) . (9) 

Mapping this quantity back to the lattice leads in fact to a sensible defini- 
tion of a Wilson loop in the NC gauge theory [25,26]. Such Wilson loops 
are complex in general [9,15]. The action ([H]) is real, however, since both 
orientations of the plaquettes are summed over, which is essential for the 
feasibility of numerical simulations. 

Simulations of gauge theories with the standard Metropolis algorithm are 
notoriously inefficient. This also holds for the TEK model. Moreover the 
usual remedy — the application of the heat-bath algorithm — cannot be 
applied straightforwardly, because the dimensionally reduced action (IHl) is 
non-linear in the link variables. However, by introducing an auxiliary matrix 
field the action can be linearised [27], so that the heat-bath algorithm works 
and the model can indeed be simulated efficientlyH This method allowed us 

■^This method has recently been extended to the 4d model where two dimensions are 
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to explore rather large systems of > 100. For the parameter sets {N,i3) 
that we investigated, we collected statistics of 1000 well thermalised and 
decorrelated configurations. 

3 The planar limit 

In the planar limit we obtain the U{N oo) lattice gauge theory on a 
commutative plane, which was solved by Gross and Witten [28]. In this limit 
they found an exact area law for the Wilson loops. In (dimensionless) lattice 
units it takes the form 

(iy(JxJ))=exp(-a(/3)/J), ^(/^) = { _ ln"(f- ^) ^ > 1/2 ' ^^^^ 

In terms of dimensional units the string tension a turns into an area, and it 
allows us therefore to identify a dimensional lattice spacing as 

a=v^. (11) 

In these units the string tension is set to 1 in the planar limit. 

Of course, we are ultimately interested in the DSL according to eq. ([7j). 
But to give it an explicit meaning we first have to identify a dimensional 
lattice spacing, i.e. we have to introduce a scale to interpret the lattice units. 
Relation ( iTTl) can be used for this purpose [15,30], provided that the values 
of A^, which are accessible to our simulations, do approximate the planar 
limit well (for the quantities of interest). Otherwise one would have to worry 
about finite N artifacts distorting the physical interpretation of our results. 

We verified this property first by checking the validity of the large 
Schwinger-Dyson equations. On the lattice, they relate Wilson loops of dif- 
ferent shapes in the planar limit An example for the corresponding contours 
is illustrated in Figure [1] (cf. first work in Ref. [22]). Indeed, we observe that 
our measurements for the two sides of this equation converge well as we 
increase A^ to a magnitude of O(IOO) at fixed /?, see Figure 

The Eguchi-Kawai equivalence to the model solved by Gross and Witten 
also implies the validity of the APD symmetry in the planar limit. Again 
we testify if this symmetry can be observed to set in (approximately) for the 

reduced; this is adequate for QED in a four dimensional space, composed of a commutative 
and an NC plane [19]. 

^This was the property that motivated the original construction of the dimcnsionally 
reduced matrix model (without twist factor) by Eguchi and Kawai [29]. 

variant of this result has been reported before in a proceeding contribution [31] and 
in a Ph.D. thesis [32]. Earlier observations in this context were given in Rcf. [30]. 
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Figure 1: An example of a set of contours, which are involved in a Schwinger- 
Dyson equation. These equations relate the vacuum expectation values of the 
corresponding Wilson loops in the planar limit. They can be derived from the 
invariance under an infinitesimal substitution of the compact link variables 
on the lattice 
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Figure 2: The convergence towards the validity of the Schwinger-Dyson equa- 
tions as N increases at fixed (3 = 0.75 (on the left) and (3 = 1 (on the right). 
The y-axis is the absolute value of the deviation between the two sides of the 
specific Schwinger-Dyson equation illustrated in FigureUl The x-axis corre- 
sponds to the lattice area of the contour C , which has a rectangular shape of 
the form I x (J + 1). These equations are well approximated as N reaches 
O(IOO). 



system sizes that we simulated. We consider four types of Wilson loops: 
we denote them as square loops, rectangular loops, stair loops and L-loops. 
In particular we considered these loops at the areas A = 4, 9, 16, 25 . . . 
in lattice units. The rectangular loops are maximally anisotropic, i.e. their 
shapes are rectangles of side lengths 1 and A^ For the shapes of stair loops 
and L-loops we refer to Figure O Note that these two loop types involve two 

^In Section 4 we will also consider rectangles with a fixed ratio between the two side 
lengths. 
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slightly different cases, depending on whether the area A is even or odd. 



Figure 3: An illustration of the different Wilson loops that we considered. 
Their contours are all polygonal, without multiple-intersections . They are 
squares, L-shapes (with legs of width 1 lattice spacing and equal length, or 
lengths deviating by 1 if the area A is even ), maximally anisotropic rectangles 
and stairs (again in two variants, depending whether A is even or odd). 

All the four types of loops (with A fixed) are related by APDs. In par- 
ticular the square and rectangles transform into each other under SL{2,R) 
on the plane. 

We now present numerical results for these Wilson loops W as we ap- 
proach the planar limit. As an example, we first consider the absolute value 
\W\ as a function of the dimensional area at a fixed lattice area of A = 36. 
Figure H] shows results for N = 75 and 155. We see that the absolute values 
for the different shapes coincide for small and for large areas, but they split 
apart at an intermediate physical area of Aa'^ = (9(10). 

However, even in this intermediate regime the differences between the 
Wilson loops at a fixed area converge to zero in the planar limit; examples 
for this behaviour are shown in Figure [5l Since f3 was kept fixed in these 
plots, both the area in lattice units and in dimensional units is constant. 

For the interpretation of our simulation results, we can therefore rely on 
the scale fill I) extracted from the planar limit. This allows us to proceed now 
to the investigation of the DSL, which describes the theory on a NC plane 
in the simultaneous UV and IR limit. 
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Figure 4: The absolute values of Wilson loops with different shapes and lattice 
area A = 36. They are plotted against the physical area 36a^, which is varied 
by using different values of (3. We show results for N = 7b on the left, and 
for N = 155 on the right. The results for different shapes coincide best at 
small and at larger areas, but they differ most around 36a^ ~ 10. 

0.01 



0.005 



-0.005 



-0.01 



-0.015 



-0.02 



1 


1 


1 1 

square - rectangle 
square - stairs 
square - L-shape 


1 

1- - * - -' 

l_. _.Q.._..J - 


S - - 


:"s:-M--_---_ 




^ -H— 
- - 


















1 


1 


A = 25, 

1 1 


0.77 _ 





0.01 


A 






0.005 






V 





o 








o 
o 


-0.005 












-0.01 








-0.015 




-0.02 



1 

H-- ._ 


1 


1 1 

square - rectangle 
square - stairs 
square - L-shape 


1 


















1 


1 


A = 36, 

1 1 


1.05 - 



60 



80 



100 
N 



120 



140 



160 



60 



80 



100 
N 



120 



140 



160 



Figure 5: An illustration of the convergence towards the APD symmetry in 
the planar limit. On the left we show the differences of {\W\) for Wilson 
loops with lattice area A = 25 at P = 0.77 , so the dimensional area amounts to 
Aa^ ~ 9.81. In the plot on the right-hand- side the parameters are A = 36 and 
(3 = 1.05, which implies almost the same area, Ac? ~ 9.79. For increasing 
N the differences between the Wilson loops with different shapes decrease 
rapidly, so we are approximating well the behaviour in the planar limit. 



4 The Double Scaling Limit 

We now employ the scale (fTTl) provided by the Gross- Witten area law ffTUl) : 
its use for the range of parameters under consideration has been justified in 
Section 3. This enables us to study the DSL (I7j) to a continuous plane of 
infinite extent and finite non-commutativity. 
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From Refs. [15] we know the following properties about the square loops: 

• The observable {W{I x J)) does indeed stabilise in the DSL. The ex- 
istence of this universality class shows in particular that the model is 
non-perturbatively renormalisable. 

• At small area, the absolute value >^ follows an area law. In 
that regime, which extends up to (a/)^<4, the phase is practically 
zero. 

• For larger areas, \{W{I x I)) \ does not decay any further, but the phase 
starts to increase linearly in the area. It obeys the simple relation 

phase = = (a/)2 ■ B , (12) 

u 

where we symbolically introduced a magnetic field B = 1/6 across the 
plane. As we mentioned in Section 2, this identification of the mag- 
netic field has been implemented in string theory and in solid state 
physics. The behaviour fll2p just corresponds to the Aharonov-Bohm 
effect, which fits perfectly into the magnetic interpretation of non- 
commutativity. Ref. [15] discovered this behaviour (unexpectedly) as 
a dynamical effect at low energy. 

To revisit the second point — the area law for small Wilson loops, where 
the phase is practically zero — we measured the Creutz ratio 

(T n-_i r (Wjlxj)) (iy((j-i)x(J-l))) ] 

'""l {W{iI-l)xJ)) {W{Ix{J-l))) \ ■ ^ 

This ratio singles out the string tension a for decays oc exp{—aA), pro- 
vided that it is equivalent for the various rectangular Wilson loops involved. 
Typical results for (nearly) square shaped Wilson loops, x(/,/), as well as 
extremely anisotropic (rectangular) Wilson loops, x(2, J), are shown in Fig- 
ure EIiI For both shapes we find a stable behaviour as we increase at fixed 
6, which suggests that our results can safely be extrapolated to the DSL. 
Deeply inside the area law regime we obtain a ~ 1. Hence in this range the 
behaviour in the DSL coincides with the planar limit. We observe, however, a 
marked deviation from it as the area approaches the transition to the regime 
of the Aharonov-Bohm type behaviour. 

At moderate area the Creutz ratios for the squares and the rectangles 
differ a little. This is a first observation hinting at shape independence for 
rectangles deeply inside the area law regime, but not beyond. 



^We actually averaged over x(2, J) and x( J, 2) in order to increase the statistics. 
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Figure 6: The Creutz ratio J) for Wilson loops with small areas, up 
to the transition at the end of the area law regime. The parameter sets are 
{N,(3) = (125,3.91), (195,6.09) and (255,7.97), wMch all corresponds to 
9 2.6, so we are approaching a DSL. For both shapes, i.e. for x(J, J) 
(squares) and for xi'^i J) (rectangles), we observe stable Creutz ratios in the 
DSL. At small area the string tension takes the same value as in the planar 
limit. At moderate area the results from square shaped and from extremely 
anisotropic rectangles begin to differ a little. 

Regarding the behaviour at large areas, which is very specific to the NC 
plane, one may wonder why the short-ranged non-commutativity has striking 
effects on the large rather than the small Wilson loops. UV/IR mixing [20] is 
apparently at work, even though the perturbative expansion of this model can 
be formulated without divergences. This suggests that UV/IR mixing occurs 
non-perturbatively, and it belongs therefore to the fundamental nature of NC 
field theory. This is in agreement with analytic [33] and numerical [32, 34] 
results for the NC model, and for 4d NC QED [19]. 

At last we proceed to a systematic study of the fate of the APD symmetry 
in the DSL. Figures [7] and [8] show results for the absolute values of Wilson loop 
expectation values, | (W) \ , with respect to the shapes illustrated in Figure 
El Figure [7] refers to a fixed non-commutativity parameter 9 = 1.63, and it 
shows results for two values of A^. For increasing A^ the dimensional volume 
is enlarged and the lattice becomes finer (because Na'^ is kept constant), so 
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we approach the DSL. The corresponding resuhs at 6 = 2.63 are presented in 
Figure [HI For the latter 6 value we show in addition the Wilson loop phases 
arg((Vr)) in Figure [9l These figures demonstrate that the Wilson loops for 

= 125 and = 155 — plotted against the dimensional area — are almost 
identical. This shows that we are indeed in the asymptotic regime of the 
DSL. We can therefore be confident that our results reveal the behaviour in 
the continuous NC plane. This confidence will be further substantiated by 
results for variants of these observables to be presented in the Section 5. 

As soon as the area exceeds the area law regime, we observe a clear dis- 
tinction between the absolute values | {W) \ for different shapes, see Figures [7] 
and [80 This distinction occurs for any pair of the contour types considered. 
We show again the rectangles of the form 1 x A (in lattice units), which be- 
come infinitesimally narrow in the DSL (in dimensional units). In addition 
we also include rectangles with a fixed ratio of 4 between the side lengths, 
which keep an invariant shape as we approach the DSL. The different expec- 
tation values for the narrow rectangles with the same area show clearly that 
the APD symmetry breaks. 

The distinction between squares and rectangles of fixed side ratio in these 
plots is less striking than the other cases. However, exactly these rectangles 
are directly relevant to explore the fate of the symmetry subgroup SL{2, R), 
hence we will focus on them specifically in Section 6. 

Next we discuss the phases arg((iy)) , which we show for 6 = 2.63 in 
Figure O As the area increases beyond the area law regime, the square 
shapes, rectangles with fixed side ratio and the stairs follow very well the 
Aharonov-Bohm type behaviour corresponding to eq. (fT2|) . which had been 
observed earlier for squares and certain rectangles [15]. The L-shape does not 
agree optimally, but its behaviour is reasonably close, in particular on the 
finer lattice which corresponds to A^ = 155. We remark that the extremely 
anisotropic rectangles considered earlier lead to strong deviations from eq. 
(fT2ll . as we are going to show in the Section 5 (Figure fT3l) . Apparently shapes 
which become extremely thin (in physical units) as we approach the DSL can 
lead to such features (although this is not the case for the stair loops). 

In any case, our results for the phases are very similar for the different A^ 
values, so they confirm that we are in an asymptotic window of the DSL. The 
phase of the Wilson loop has a much stronger trend towards (at least partial) 
APD symmetry than the absolute value. As far as we could check, the phase 
(alone) is well compatible with this symmetry for fixed shapes in dimensional 
units, with a finite extent in each direction. But of course the differences in 

^We note that the total area of the system amounts to V = ttNO, hence for our 
parameters V ^ 9 is granted. 
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Figure 7: The absolute values of various Wilson loops at a fixed non- 
commutativity parameter 9 = 1.63. On top we show results at N = 125 
and below at N = 155. As we increase N the dimensional volume grows and 
the lattice spacing shrinks, so that we approach simultaneously the limits to 
the continuum and to infinite volume (UV and IR limit). The striking sim- 
ilarity of these plots confirms that DSL convergence is reached. The results 
show a clear shape dependence beyond the area law regime, and therefore the 
breaking of APD symmetry. 
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Figure 8: The analogous plots to Figure\7\ but now at 6 = 2.63. Again the 
results clearly confirm a shape dependence beyond the area law regime, which 
is stable as we increase N towards the DSL. 



\{W)\ are sufficient to discard APD symmetry as a basic property of this 
theory. 
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Figure 9: The phases of various Wilson loops at 9 = 2.63 for N = 125 (on 
top) and N = 155 (below). The Aharonov-Bohm type behaviour in eq. fT^j 
at large area is very well confirmed for the shapes shown here. They include 
results for rectangular loops with a fixed ratio of A and a ratio of 11 between 
the side lengths, so that its shape is invariant on the way to the DSL. The 
L-shape, on the other hand, becomes infinitesimally thin in this limit. Its 
phase follows eq. (W) with some fluctuations. 



15 



5 APD symmetry breaking for alternative ob- 
servables 



Our results for the standard observables | (W) \ and arg((Vr)) were presented 
in Section 4 and illustrated the shape dependence of Wilson loops with the 
same area (in dimensional units). We gave evidence for this effect to persist 
in the DSL. 

In this section we present the corresponding results for {\W\) and {aig{W)), 
which are somewhat different observables. The quantities of Section 4 are 
tractable in perturbation theory, and these are therefore the observables that 
have been addressed in Refs. [8-14]. Moreover they were measured in the 
previous numerical study in Ref. [15]. On the other hand, the quantities of 
this section are generally not considered in analytic work. In the appendix we 
comment on the prospects of the perturbative treatment of However, 
{\W\) and {a.Tg(W)) can be handled numerically without specific problems, 
and they also represent valid physical observables]^ For instance, in lattice 
gauge theory it is also usual to measure {\P\) {P being the Polyakov loop) — 
its magnitude serves as a criterion to distinguish the phases of confinement 
and deconfinement. 

We add these observables here in order to supply further strength to our 
observation of APD symmetry breaking. The new observables are suitable 
for this purpose, in particular because does not become as tiny as 

\{W)\ at moderate and large area, which leads to smaller relative errors. 

We first present the results for {\W\) for two values of at 6^ = 1.63 
(Figure [TO!) and at 6* = 2.63 (Figure [TTI). in analogy to Figures [7] and [H The 
new observables are practically identical to those of Section 4 at small area, 
where the phase is tiny throughout the Monte Carlo history. They differ, 
however, as the area grows beyond this regime. There is significantly 

larger, which enables a better distinction (without overlapping error bars). 
But the qualitative behaviour is very similar to Section 4. 

Figure [12] shows the phases {aTg{W)) at the same parameters as in Figure 
[TUl It can be compared to the phase a.Tg{{W)) presented before in Figure 
[9] (for 6 = 2.63). Also this phase follows closely the Aharonov-Bohm type 
behaviour of eq. ( fT2|) for the squares, the stairs loops and for the rectangle 
with a fixed side ratios. (For the squares this behaviour was also observed 
in 4d NC U{1) gauge theory [19].) Again the phase of the L-shape loops 
deviates. It fluctuates strongly in the Monte Carlo history all over the interval 
(— TT, n] , keeping its mean values close to zero. 

^Actually wc considered the observables of this section ah'eady in Section 3 when we 
discussed the vahdity of the scale identified from the planar hmit (Figures S] and [S]). 
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Figure 10: The expectation values {\W\) for various Wilson loops at a fixed 
non-commutativity parameter 9 = 1.63. On top we show results at N = 125 
and below at N = 155. The striking similarity of these plots confirms that 
DSL convergence is reached. The results reveal a clear shape dependence 
beyond the area law regime also for this observable. 

At last we consider the phases of the rectangles of the form 1 x A, i.e. 
of an infinitesimally narrow shape in the DSL. We mentioned in Section 4 
that these phases do not follow eq. (fT2l) . In Figure [13] we add the phases for 
these rectangles at ^ = 2.62. We show arg((W^)) on the left, and {aig{W)) 
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Figure 11: The analogous plots to Figure [70| but now atO = 2.63. The results 
affirm once more a shape dependence beyond the area law regime, which is 
stable as we increase N . In particular the differences at fixed area do not 
shrink as we increase N towards the DSL, in contrast to the planar limit 
behaviour in Figure 

on the right. They are very similar, and — more importantly — they are in 
excellent agreement for the different values. This confirms once more that 
we access the double scaling window with the parameters used. 
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Figure 12: The phases (arg(H^)) of various Wilson loops at 9 = 1.63 for 
N = 125 (on top) and N = 155 (below). The Aharonov-Bohm type behaviour 
in eq. [W) at large area holds well also for this phase for the stair loops, the 
squares and the rectangular loops with a fixed side ratios of 4 and of 11. On 
the other hand, the average phase of the L-shape loops remains very small 
even at large areas. 



This section presented additional APD symmetry breaking results, now 
for {\W\) and (arg(iy)). The values and physical meanings of these observ- 
ables are different from those in Section 4, but qualitatively the results are in 
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We show the phase of (W) (on the left) and the expectation value of arg(iy) 
(on the right). Their values are very similar, but these phases are not related 
to eq. ^W) . The essential observation here is that the phases practically 
coincide for the two sets of parameters, which confirms again that we do see 
the window that extrapolates to the DSL. 



full agreement. The observables added here provide further strength to our 
investigation of APD symmetry breaking — see also the theoretical reasoning 
in the Appendix. 



6 The 5L(2, R) symmetry breaking 

While the breaking of the full APD symmetry has already been demonstrated 
extensively, the specific case of the SL{2, R) symmetry subgroup may seem 
less obvious from the Figures shown so far. Hence this section focuses on 
rectangles only to illustrate in particular the breaking of this subgroup. 

Figures [m and [13 are dedicated to the decay of the absolute values 
and for the rectangles. We show in both figures the behaviour for fixed 

6, where we include the data from = 125 and A^ = 155 in the same plot. 
Here we also add results for the decay of rectangles with a fixed side ratio of 
11, which are helpful to demonstrate the SL(2,R) breaking more clearly. 

For the rectangles with fixed shapes these observables tend to oscillate as 
we vary the dimensional area beyond the area law regime. This observation 
is well compatible with the data from both A^ values, hence it seems to 
characterise the DSL. These oscillations have similar mean values, but the 
amplitude is significantly larger for the ratio 11 between the side length when 
compared to ratio 4 or 1 (squares). 

We conclude that the SL{2, R) symmetry is indeed more viable to some 
approximation than the rest of the APD symmetry group, but its breaking 
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is nevertheless manifest on the non-perturbative leveL 
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Figure 14: The decay of the absolute values \{W)\ (on top) and of {\W\) 
(below) for rectangular loops at 6 = 1.63. Beyond the area law regime, the 
rectangles with finite sides in the DSL (i.e. with a fixed side ratio) tend to 
oscillate around ~ 0.01 for \{W)\ and a somewhat larger value for 
The amplitudes, however, depend clearly on the side ratio, which shows the 
SL(2, R) symmetry breaking. This is most evident for the rectangles with 
side ratio 11 (their data are connected by a line to guide the eye). 



21 



10 



0.1 



0.01 



0.001 



1e-04 



e = 2.63 



N=125: squares r- -+-- -; 
rectangles (1 : 4) 
rectangles (1 : 1 1) a 
rectangles (1 x A) ^ q 
N=155: squares ^ x ; 
rectangles (1 : 4) 
rectangles (1 : 1 1) i- 
rectangles (1 x A) ■ 







4 






M + 


^ • 


\ 

•CD 


■- 


X ^_ 


r'" 


tf : 
m'' i 

ti ; 

1 ; 


25 


30 



dimensional area 



A 

5 

V 



0.1 



0.01 



N=1 25: rectangles (1 X A) --a---: 

squares ^ + : 
9 = 2-63 rectangles (1 : 4) 

rectangles (1:11)---' 

N = 155: rectangles (1 x A) ^ ■ 

squares ^ x : 

rectangles (1 : 4) i- -i 

rectangles (1 : 1 1) ^- ■» 



i 



* * x " ■ / - 



10 15 20 25 30 35 40 
dimensional area 



Figure 15: The analogous plots to Figure but now at 9 = 2.63. The 
observed decays of the rectangular Wilson loops is very similar for both non- 
commutativity parameters 9 that we consider. Again the SL{2, R) symmetry 
breaking can been seen best from the data for rectangles with side ratio 11 
(which are connected to guide the eye). 

7 Conclusions 

While APD symmetry holds in ordinary two-dimensional Yang-Mills theory, 
it turns out to be broken by non-commutativity. Perturbation theory can 
reveal this effect, depending on the order considered. Since this is still re- 
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strictive, to a finite order the symmetry may be protected in part [11], with 
a residual subgroup showing up in other approaches as well [13]. 

However, on the non-perturbative level and by considering finite 9, the 
quantum effects unfold their full power and destroy the APD symmetry, 
including the subgroup SL{2,R). Here we presented explicit results for 
the APD symmetry breaking. We considered four types of Wilson loops 
with polygonal contours, without multiple-intersections. At the same area, 
numerical simulations revealed the shape dependence of the Wilson loop 
expectation values on the lattice. Our results were obtained for constant 
non-commutativity parameters 6 at different values of — for increasing 
A^ the volume becomes larger (in dimensional units) and the lattice finer. 
The results for the Wilson loops remain very stable, hence they allow for 
a reliable extrapolation to a continuous plane of infinite extent at fixed 
non-commutativity (Double Scaling Limit). This limit reveals the break- 
ing of APD symmetry — including the SL{2, R) subgroup — on the non- 
perturbative level. Furthermore our results extend the loss of APD symmetry 
to the discretised model in a finite volume. 

To be more precise, the APD symmetry holds very well for relatively small 
loop areas (in dimensional units). For a fixed lattice area this corresponds to 
the weak coupling regime, and hence the agreement with the leading order 
perturbative results appears consistent. However, as the dimensional area 
increases, this agreement — and therefore the APD symmetry — collapses: 
it partially persists in the phase of the Wilson loops, but not in their abso- 
lute values. We also observed that the SL(2, R) subgroup has is a relatively 
high viability as an approximate symmetry (which is consistent with pertur- 
bative results), but its breaking becomes manifest for strongly anisotropic 
transformations in this subgroup. 

We conclude that NC gauge theory has a rich structure, even in = 2, 
far beyond its commutative counterpart. This conclusion rises hope for ex- 
citing effects. On the other hand, it also means that an analytical solution 
is unlikely. Therefore numerical results are of great importance, as it is the 
case in commutative Yang-Mills theories in four dimensions. 
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A The perturbative treatment of {\W\) 

In this appendix we comment on the perturbative treatment of the observable 
along the lines of the formalism constructed in Refs. [8,9]. We refer 
the reader to those works for the details and the relevant literature, while 
here we just recall the basic ingredients needed for our argument. 

A convenient gauge choice for perturbative Wilson loop calculations is 
the light-cone gauge A_ = 0. Faddeev-Popov ghosts are known to decouple 
also in the NC case, and the two-dimensional Lagrangian in this gauge looks 
indeed free. The non-trivial information on the dynamics is encoded in the 
singular behaviour of the propagator 

D^^ = i[k-J]. (14) 

Two prescriptions for the pole are important, the Cauchy principal value 
method by 't Hooft, and the prescription due to Wu, Mandelstam and Leib- 
brandt [35], which employs the propagator 

D++ = i[k- + iek+]-\ (15) 

In the commutative case, the latter prescription is genuinely perturbative, 
while the former is able to produce results for the Wilson loops, which take 
into account non-perturbative contributions, and it yields exponentiation 
(area law). But in the NC case, the former prescription behaves very wildly 
and has been soon abandoned, while the latter gives sensible results at all the 
orders considered. Also, the rotation to the Euclidean version of the theory 
turned out to be useful. 

When examining the perturbative series of the Wilson loop (W), one 
notices that the NC phase factor intermingles non-trivially, in the non-planar 
diagrams, with the propagators in the momentum and contour integrations. 
The results are themselves expansions, typically in 1/9; the Wilson loop is 
analytic around 6 = oo. 

The perturbative series for the observable {\W\) looks a priori quite dif- 
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ferent from the one of The relevant expansion in the U{1) case reads 



/ I °° 

\ n=0 ^ J s„-i 

j d?x + ^ . . . * + e(s„)) I ^ . (16) 

Clearly the presence of the modulus inside the quantum average complicates 
the task of rearranging the expansion in terms of the Green functions of the 
theory, for one has first to write it schematically as \/ (^ --OIX] •••)* then 
perform a Taylor series in g. The fact that the sums begin with a constant, 
followed by terms with an increasing number of fields, allows to formally 
treat the series as analytic around g = 0. At the end, one computes the 
vacuum expectation value, based on the Green functions 

(0|TA+(x + e(si))*...*A+(x + e(s2n))|0) . (17) 

Each term can in principle be treated similarly to the standard case: one 
is left with multiple integrals of propagators and Moyal phases, in momen- 
tum space and on domains which are simplexes in the contour parameters. 
These contributions can be visualised by attaching the propagator lines along 
the contour itself: different simplexes correspond then to different crossing 
patterns of the propagators. 

In the simplest case of the 0{g^) computation, the momentum integrals 
can be carried out exactly using complex plane techniques and identities of 
Bessel functions, and the remaining contour integration is then performed 
numerically. For the new observable, care must be taken in computing some 
of the integrals, because one encounters terms like 

(^J dsii^{si) j (fx exp[ip{x + ^{si))]A+{p) 

which appear already at the order 0{g'^). This order is unaffected by the 
Moyal phase, but still the overall integration over the space-time base-point 
X forces the appearance of a correlator {A^{0)A_^.{0)) at zero momentum 
which is infrared divergent. On the other hand, the integration over the 
closed contour vanishes because it reduces to dsiC,-{si) = 0. Therefore, 
one has to fix the order of integrations, such that the contour integration is 
performed first. Still, it is not clear if this can be made consistent at any 
order, or if other potentially divergent terms will appear due to the nesting 
of contour integrals, or, most important, due to interference with the NC 
phase factor. We leave this interesting analysis for future developments. 
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Here we would like to remark that, provided one can consistently treat these 
ambiguities, then the arguments about order-by-order APD invariance at 
6 = seem to apply. This implies that the observed shape- dependence of 
{\W\) in the NC plane appears as a valid argument for concluding that APD 
symmetry breaks down when going from 6' = to a finite 6 value. 
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